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Abstract: Intense human activities in arid areas have great impacts on groundwater hydrochemical cycling 
by causing groundwater salinization. The spatiotemporal distributions of groundwater hydrochemistry are 
crucial for studying groundwater salt migration, and also vital to understand hydrological and 
hydrogeochemical processes of groundwater in arid inland oasis areas. However, due to constraints posed 
by the paucity of observation data and intense human activities, these processes are not well known in the 
dried-up river oases of arid areas. Here, we examined spatiotemporal variations and evolution of 
groundwater hydrochemistry using data from 199 water samples collected in the Wei-Ku Oasis, a typical 
arid inland oasis in Tarim Basin of Central Asia. As findings, groundwater hydrochemistry showed a 
spatiotemporal dynamic, while its spatial distribution was complex. TDS and 6!8O of river water in the 
upstream increased from west to east, whereas ion concentrations of shallow groundwater increased from 
northwest to southeast. Higher TDS was detected in spring for shallow groundwater and in summer for 
middle groundwater. Pronounced spatiotemporal heterogeneity demonstrated the impacts of geogenic, 
climatic, and anthropogenic conditions. For that, hydrochemical evolution of phreatic groundwater was 
primarily controlled by rock dominance and evaporation-crystallization process. Agricultural irrigation and 
drainage, land cover change, and groundwater extraction reshaped the spatiotemporal patterns of 
groundwater hydrochemistry. Groundwater overexploitation altered the leaking direction between the 
aquifers, causing the interaction between saltwater and freshwater and the deterioration of groundwater 
environment. These findings could provide an insight into groundwater salt migration under human 
activities, and hence be significant in groundwater quality management in arid inland oasis areas. 
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1 Introduction 


Groundwater plays an important role in the socio-economic development and eco- 
environmental protection in arid areas due to scarce precipitation and deficient surface water 
(Lezzaik et al., 
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2018; Erler et al., 2019; Jasechko et al., 2021), especially in arid inland oasis areas (Wang et al., 
2021). In the last decades, groundwater depletion and salinization in arid areas have led to a series 
of hydrological and ecological havocs owing to climatic variations and human activities that 
challenge the local water and ecological security (Zhang et al., 2014; Asoka et al., 2017; de Graaf 
et al., 2019; Chen et al., 2020; Kulmatov et al., 2020; Zhao et al., 2021). Still, groundwater 
salinization can affect public health, soil fertility, vegetation growth, and renewable groundwater 
resources (Tweed et al., 2011; Fuchs et al., 2019; Jia et al., 2020). Groundwater hydrochemistry is 
very important in biogeochemical cycling as it may influence the freshwater utilization and 
ecological functioning of surface water and groundwater in arid areas (Cary et al., 2015; Obeidat 
et al., 2021). Additionally, the spatial distributions of dissolved ion concentrations in groundwater 
could provide insight into the origin, migration pathways, and evolution of groundwater salt (Liu 
et al., 2018; Wang et al., 2018), which are helpful to understand the hydrological and hydro- 
chemical cycling in arid areas (Fuchs et al., 2019; Wu et al., 2020; Shen et al., 2021). Therefore, 
there is a need to improve our understanding of the hydrochemical characteristics and evolution 
of groundwater influenced by intense human activities in arid inland oasis areas. 

The Tarim River Basin, the largest inland river basin in Xinjiang Uygur Autonomous Region of 
China, is located in the arid inland area of Central Asia and is also the core area of constructing 
the "Silk Road Economic Belt" (Chen et al., 2019; Wang et al., 2021). In the last 30 years, with a 
continuous intensification of farmland expansion and population growth, groundwater 
overexploitation and the increased surface water diversion from the river for irrigation purposes 
have caused the river to dry up, phreatic decline, groundwater environment deterioration, 
ecological degradation, and desertification in the oasis-desert regions of the Tarim Basin (Wang et 
al., 2013; Xiao et al., 2015; Wang et al., 2020). Accordingly, the deterioration of the groundwater 
environment in arid oases has influenced hydrogeochemical and ecological processes, plant 
communities, and soil salinization (Liu et al., 2018; Jia et al., 2020; Noshadi et al., 2020). Thus, 
the ecological restoration of arid inland oasis areas should be based on the improvement of the 
groundwater environment, by focusing on the spatiotemporal characteristics and evolution 
mechanisms of groundwater hydrochemistry (Wang et al., 2013). Nevertheless, due to observation 
data scarcity and the complex interaction between surface water and groundwater following 
intense human activities, the hydrochemical evolution and migration mechanism of groundwater 
remain poorly understood in the Tarim River Basin. All these may affect the implementation of 
sustainable groundwater resources management and ecological restoration (Chen et al., 2019; 
Zeng et al., 2020). 

Based on the hydrochemical and isotopic techniques, remote sensing inversion, and 
hydrological models, previous studies in arid inland oasis areas mainly focused on the origin and 
age (Ma et al., 2013), spatial distribution and temporal variation (Gates et al., 2008; Hasan et al., 
2011; Pant et al., 2018), and evolution mechanism and migration (Wang et al., 2013; Huang et al., 
2017) of hydrochemical components for phreatic groundwater and confined water. As well, the 
fresh-salt water interaction between phreatic groundwater and confined water has been analyzed 
(Liu et al., 2018; Zeng et al., 2020). Last but not least, the impacts of climate change and human 
activities on groundwater hydrochemistry in arid inland areas have also been studied (Xiao et al., 
2015; Wang et al., 2021). In addition, several mechanisms could control the groundwater 
hydrochemistry, including the flushing, evaporation-condensation, dissolution and leaching, 
cation exchange-adsorption, and mixing (Xiao et al., 2015; Wang et al., 2018; Li et al., 2019; 
Zeng et al., 2020). Groundwater hydrochemistry in arid inland areas may be affected by the 
climatic and hydrogeological conditions, groundwater recharge, land cover types, groundwater 
exploitation, and agricultural activities among others (Scanlon et al., 2006; Kaur et al., 2019; 
Riley et al., 2019; Jia et al., 2020). Yet, the spatiotemporal heterogeneity and dominant factors of 
groundwater hydrochemistry have not been thoroughly investigated in the dried-up river oasis of 
Tarim Basin, particularly the impact of intense human activities. 

Here, we analyzed the spatiotemporal variations and evolution of groundwater hydrochemistry 
using data of 199 water samples collected from August 2018 to August 2019 in the Wei-Ku Oasis, 
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a typical dried-up river oasis of the Tarim Basin in Central Asia. The specific objectives of this 
study are to: (1) examine the spatiotemporal variations of groundwater hydrochemistry; (2) 
analyze the dominant factors influencing groundwater hydrochemistry evolution; and (3) explore 
groundwater salt migration under the influence of human activities. The results are expected not 
only to enhance our understanding of groundwater salt migration under intense human activities 
but also to help improve groundwater management in the Tarim Basin, hence the prevention of 
groundwater environment deterioration. 


2 Materials and methods 


2.1 Study area 


The Tarim River Basin (73°24’-93°39’E, 34°48'-43°21'N), surrounded by the Tianshan 
Mountains and the Kunlun Mountains in Xinjiang Uygur Autonomous Region, covers an area of 
1.02x10° km? and is the largest inland river basin of Northwest China (Fig. 1). The Wei-Ku Oasis 
(82°05'—83°28'E, 40°54'-41°51'N), located in the middle and lower reaches of the Weigan-Kuqa 
River Basin in the northern Tarim Basin, covers an area of 7.10x10° km’, with the elevation 
ranging from 945 to 1147 m above sea level. It mainly includes the Weigan River and Kuqa 
River; the Weigan River is composed of the Kizil River, Muzhate River, Televichuk River, 
Kapusilang River, and Karasu River. The study region is characterized by a temperate arid 
continental climate. The annual mean air temperature and mean annual precipitation values are 
11.4°C and 74.6 mm, respectively (data from the Chinese National Meteorological Center; 
http://data. cma.cn/). The highest air temperature is generally recorded in July or August, and the 
lowest is in January (Fig. 2). Precipitation mainly occurs from May to September, accounting for 
more than 76% of the total annual precipitation, and peak rainfall occurs in June or July. The 
annual runoff flowing into the Wei-Ku Oasis is about 31.2x10° m? (more than 90% of the runoff 
using for irrigation); it has three peaks due to reservoir regulation, occurring in March, July, and 
November. The groundwater level depth of the oasis ranges from 2 to 6 m, and groundwater level 
generally declines from northwest to southeast (Wang et al., 2021). 

The Wei-Ku Oasis is a plain region with unconsolidated Quaternary sediments, which 
comprises the interstratified pebbly sandstone, fine sandstone, sandy gravel, silty fine sandstone, 
sandstone, and fine siltstone, forming the porous aquifers and aquitards (Fig. 1c) (Wang et al., 
2021). Within the oasis, from the northern piedmont plain to the southern alluvial plain, the 
particle size of the rock ranges from coarse to fine, the aquifer system gradually transits from 
single structure to multi-layer structure, and the groundwater level depth declines gradually. The 
aquifer sediments are pebbly sandstone and sandy gravel in the piedmont plain. The pebbly 
sandstone transits to silty fine sandstone in the middle alluvial-diluvial plain and fine sandstone in 
the alluvial plain (Fig. 1c; Wang et al., 2021). Furthermore, the dominant natural plant species are 
Sophora alopecuroides, Tamarix spp., Populus euphratica, Karelinia caspica, Phragmites 
communis, Halostachys caspica, Glycyrrhiza sp., and Alhagi pseudalhagi. The land cover types 
within the Wei-Ku Oasis were mainly dominated by cropland (62.5% in this oasis), grassland 
(17.9%), and bare land (13.0%) in 2018. The cereal crops within the oasis mainly included wheat 
and corn, and dominant economic crops were cotton, sugar beet, and rapeseed (Chen et al., 2019). 


2.2 Data and methods 


Surface water and groundwater samples were collected from August 2018 to August 2019 along 
the Weigan River and Kuga River, as well as within the Wei-Ku Oasis (Fig. 1b), which falls in 
three seasons (spring, summer, and autumn). During the sampling period, precipitation samples 
were also collected after each day of precipitation (rainfall and snowfall) from two sites in the 
upper reaches of the Weigan-Kuga River Basin (Fig. 1b). River water samples were collected 
from the upstream of Weigan River and Kuga River, along with 6 reservoir water samples that 
were obtained from the oasis. Additionally, groundwater was sampled from pumped wells 
(domestic wells, irrigation wells, and industrial wells) and regularly used boreholes within the 
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Wei-Ku Oasis. We grouped groundwater samples into three classes based on hydrogeological 
conditions (Wang et al., 2013), namely shallow groundwater (well depth lower than 20 m), middle 
groundwater (well depth ranging from 20 to 100 m), and deep groundwater (well depth greater 
than 100 m). 
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Fig. 1 Overview of the Tarim Basin and location of the Wei-Ku Oasis (a), distribution of water sampling sites in 
the Weigan-Kugqa River Basin (b), and hydrogeologic conditions of the Wei-Ku Oasis (c; modified from Wang et 
al. (2021)). 
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Fig. 2 Intra-annual variations of monthly average air temperature and monthly precipitation in the Wei-Ku Oasis 
(a) and monthly runoff in the upstream of the Weigan River (b) 


A total of 199 water samples were collected for hydrochemical and stable isotope analyses, 
including 23 precipitation samples, 101 surface water samples (95 river water samples and 6 
reservoir water samples), and 75 groundwater samples (53 shallow groundwater samples, 21 
middle groundwater samples, and 1 deep groundwater samples). All samples were first filtered 
through a 0.22-um filter, then immediately kept in 100 mL polyethylene bottles, wrapped with 
parafilm for avoiding evaporation fractionation, and lately stored in the refrigerator at 4°C until 
laboratory analyses. 

All the hydrochemical components and stable isotopes of the collected samples were measured 
at the State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and 
Geography, Chinese Academy of Sciences in Xinjiang of China. Major cations and anions 
(chloride, sulfate, sodium, potassium, magnesium, and calcium) of samples were determined by 
an ion chromatograph (Dionex ICS-5000, Thermo Fisher Scientific, Massachusetts, USA), while 
HCO; and CO,” were analyzed using the titration method. The analytical precision of major ions 
was below 1%, with the detection limit of 0.1 mg/L. The ionic charge balance errors for all water 
samples were within +10%. Electrical conductivity (EC) and temperature were measured in the 
fields using a conductivity meter (YSI ProPlus, Yellow Springs Instrument Inc., Ohio, USA), with 
an EC measurement precision of 0.01 uS/cm (Wang et al., 2020). In addition, the 3'°O and 8H of 
water samples were measured using a liquid water isotope analyzer (LGR DLT-100, Los Gatos 
Research Inc., California, USA) and reported relative to the Vienna Standard Mean Ocean Water 
(VSMOW) in delta per million (8, %o). Measurement precisions for 5°H and 8'°O were +0.8%o 
and +0.1%o, respectively. 

In this study, meteorological data (air temperature and precipitation) were obtained from the 
Chinese National Meteorological Information Centre (http://cdc.cma.gov.cn), and runoff data for 
the outlets of the mountainous tributaries were collected from the Weigan River Basin Authority, 
Xinjiang, China. Besides, data statistical analyses were conducted using the SPSS software 
package (version 22.0), ArcGIS (version 10.5), MATLAB (version R2018a), and Origin (version 
2017). The inverse distance weighting method was used to interpolate the hydrogeochemical and 
meteorological data into 250 m grid data. The Piper (1944) and Gibbs (1970) plots were used to 
further elaborate the findings. 


3 Results 


3.1 General characterization of hydrochemistry 


The stable isotope values (8°H and 5'*O) for precipitation, surface water, and groundwater are 
shown in Figure 3. The 5°H and 8'%O values of different water components presented a large 
range in our study area and differed significantly among various water components. In general, 
the stable isotope values become more depleted successively from reservoir water to 
precipitation, river water, and groundwater, which may be attributed to the difference in 
evaporation and water sources (Wang et al., 2020). Both the 5O and 8°H values of precipitation 
exhibited a relatively wide range, varying from —15.8%bo to 0.1%o and from —116.3%o0 to 11.5%bo, 
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respectively. The mean isotopes of precipitation (8'SO= —7.6%o and &° H= —48.0%o) in our study 
area were higher than those in the Urumqi station in northern Xinjiang (580 = —10.7%o and 5°H= 
—72.8%0, which were from the Global Network of Isotopes in Precipitation (GNIP)), but lower 
than those in the Hotan station in southern Xinjiang (8'*O=—5.5%o and 8° H= —32.6%o, which were 
from the GNIP), indicating the complex water vapor sources and local recycling processes (Chen 
et al., 2019). 

Surface water data were located near the global meteoric water line (GMWL; &°H=85'*O+10; 
Craig, 1961), while most groundwater samples were located above the GMWL (Fig. 3). The 
stable isotope values of river water (mean 3°H= —53.6%o and mean 8'8O= —8.9%o) during the 
sampling period were lower than the values of precipitation, maybe due to the recharge from 
glacier/snow meltwater and groundwater in the upstream (Chen et al., 2019). The isotope values 
of reservoir water ranged from —9.3%o to —5.2%o for 5°O (mean o= —6.5%o0) and from —57.9%o 
to —35.3%o for 8°H (mean H= —42.6%o). These were higher than the values of river water, due to 
strong evaporation fractionation in such arid region (Wang et al., 2013). Furthermore, the mean 
values of 8'°O and 8H for shallow groundwater (from —12.2%o to —7.5%o and from —83.9%o to 
—49.7%o, respectively) were similar to the mean values of middle groundwater (5'%O ranging from 
—10.3%o to —8.4%o and &°H ranging from —73.9%o to —59.3%o) and deep groundwater (5'°O= 
~9.5%o and &°H= —61.2%o) in the Wei-Ku Oasis. 
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Fig. 3 Stable isotopic compositions for precipitation, surface water, and groundwater samples in the Wei-Ku 
Oasis. Precipitation isotope data were from GNIP (Global Network of Isotopes in Precipitation) for the Hotan 
(HT) and Urumqi stations. WK, the Weigan-Kuqa River Basin. GMWL, global meteoric water line (Craig, 1961). 


Table 1 presents the hydrochemical results of precipitation, surface water, and groundwater. 
Generally, main dissolved ion concentrations in groundwater were significantly higher than those 
in surface water, whereas ion concentrations in precipitation were lower compared to those in 
surface water (CI, S04”, HCO;, Na, K’, Mg”, and Ca”; Table 1), mainly due to strong 
evaporative concentration, mineral dissolution, and leaching (Jia et al., 2020). Precipitation 
samples had low ion concentrations and salinity (Cl =4.5 mg/L and TDS (total dissolved 
solids)=91.7 mg/L), with a hydrochemical type of HCO3-Ca (Table 1; Fig. 4a). River water 
samples were slightly alkaline (mean pH=7.95) and were freshwater (TDS=379.4 mg/L). The CI 
concentration of river water samples ranged from 1.6 to 146.5 mg/L, and the EC varied from 208 
to 1206 us/cm, with SO4,-HCO3-Ca chemistry. Compared to river water, reservoir water in the 
oasis had slightly higher ion concentrations and salinity (EC ranging from 665 to 1361 us/cm and 
Cl ranging from 67.7 to 173.3 mg/L), with identical chemistry (HCO3-SO,-Ca), which indicates 
that reservoir water was directly from the river but undergone strong evaporation. 

Samples of middle groundwater were intermediate between the samples of shallow groundwater 
and deep groundwater (Fig. 4b), which provided further evidence for the close water relationship 
between various aquifers in this oasis. The pH of groundwater in the oasis generally ranged from 
7.48 to 8.42, indicating a slightly alkaline character, which resulted from the extremely arid 
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condition and slow groundwater flow in our study region (Wang et al., 2013). Groundwater 
samples exhibited a high degree of hydrochemical variability, and the TDS ranged from 314.0 to 
24,156.4 mg/L (varying from freshwater to saline water), suggesting the great variation of 
groundwater quality. Vertically, highly saline groundwater was observed in the shallow aquifer, 
whereas lowly saline groundwater was mostly found in the deep aquifer (Table 1). Most shallow 
groundwater had relative high salinity and ion concentrations (mean TDS=7560.3 mg/L and mean 
Cl =2875.2 mg/L), with a hydrochemical type of Cl-SO,-Na-Mg. Middle groundwater in the oasis 
had a similar hydrochemical type (Cl-SO,4-Na-Mg-Ca) compared to shallow groundwater, with an 
average TDS of 1279.5 mg/L, implying the direct recharge from the overlying shallow aquifer 
(Zeng et al., 2020). Besides, deep groundwater was dominated by Na and Ca cations and HCO; 
and SO, anions, with low ion concentrations and salinity (TDS=314.0 mg/L), and 
HCO3-SO,-Cl-Na-Ca chemistry (Table 1; Fig. 4b). It was obvious that there was a large difference 
in ion concentrations between the shallow, middle, and deep groundwater samples in the Wei-Ku 
Oasis, revealing that saline shallow groundwater could leak downward to confined water and 
undergo a mixing process between salt and fresh water (Liu et al., 2018). 


Table 1 Hydrochemical characteristics for the surface water and groundwater samples obtained in the 
Weigan-Kuga River Basin 


n Statistic pH EC cI SO," HCO; Na” kt Mg” Ca?" Hydrochemical 
(us/cm) (mg/L) type 
Precipitation 
Max. 12.6 34.7 136.2 10.1 1.9 2.6 51.2 
23 Min. 0.6 1.1 4.4 0.0 0.1 0.3 3.6 HCO;-Ca 
Mean 4.5 9.2 40.7 2.8 0.7 1.1 21.1 
River water 
Max. 8.25 1206 146.5 269.6 372.5 92.6 9.2 38.6 138.7 
95 Min. 7.48 208 1.6 23.7 75.2 22 1.0 6.0 44.5 SO4-HCO3-Ca 
Mean 7.95 609 46.2 122.9 145.9 33.0 44 20.7 79.2 
Reservoir water 
Max. 7.81 1361 173.3 259.7 244.8 110.4 10.0 46.6 113.7 
6 Min. 7.55 665 67.7 118.3 61.5 61.4 6.9 23.9 48.5 HCO3-SO,-Ca 
Mean 7.73 935 112.4 175.0 137.8 84.5 8.3 31.3 80.3 


Shallow unconfined groundwater (well depth<20 m) 
Max. 8.14 37,400 11,703.6 5502.2 754.7 6732.7 88.8 1085.6 1377.2 

53 Min. 7.52 1134 97.0 183.7 79.0 107.9 0.0 61.2 12.2 Cl-S0O,-Na-Mg 
Mean 7.83 11,477 2875.2 1674.4 358.6 2018.7 23.4 332.9 456.3 

Middle unconfined/confined groundwater (well depth in the range of 20-100 m) 
Max. 8.42 4520 686.8 1164.1 543.6 452.2 15.8 200.0 291.7 

21 Min. 7.48 1012 119.1 120.9 110.5 75.2 0.3 13:7 27.1 Cl-SO,-Na-Mg-Ca 
Mean 7.84 2033 272.4 434.4 261.4 197.3 8.3 83.5 152.8 

Deep confined groundwater (well depth>100 m) 

1 8.07 572 62.4 87.0 116.9 50.1 3.5 16.3 36.3. HCO;-SO,-Cl-Na-Ca 


Note: n, number of water samples; EC, electrical conductivity; Max., maximum value; Min., minimum value. 


3.2 Spatial variations of hydrochemistry 


Generally, the salinity of river water and shallow groundwater increased along the flow paths 
from the recharge to discharge area in our study region (Figs. 5 and 6), connoting the occurrence 
of evaporation (Wang et al., 2013; Liu et al., 2018). Figure 5 showed the spatial distributions of 
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Fig. 4 Piper plots for precipitation and surface water samples (a) as well as groundwater samples (b) in the 
Wei-Ku Oasis 
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Fig. 5 Spatial distributions of river water hydrochemistry, air temperature and precipitation in the upstream of 
Weigan-Kuqa River Basin. (a), TDS of river water; (b), 8'50 of river water; (c), mean annual air temperature 


(2018-2019); (d), mean annual precipitation (2018-2019). TDS, total dissolved solids; Tem., air temperature; 
Pre., precipitation. 


TDS and O for river water in the upstream of the Weigan-Kuqa River Basin. The TDS and 
5'°O values of river water presented obvious spatial heterogeneity. The 5'°O of river water in the 
upstream ranged from —10.9%o to —6.4%o; the values were much lower in the western part than in 
the eastern part, most probably due to much larger snow area in the west, resulting in relatively 
low isotopes of river water recharge from snow meltwater (with low stable isotope values) (Wang 
et al., 2020). Similarly, the TDS values of river water ranged from 133.7 to 680.9 mg/L and the 
values were lower in the west (Muzhate River) than in the east (Kuqa River). Overall, the TDS 
concentration of river water in the upstream increased from west to east. 

As shown in Figures 6 and 7, hydrochemical components of shallow and middle groundwater 
across the Wei-Ku Oasis varied among the sampling sites, exhibiting spatial heterogeneity, 
probably due to lithology and underlying surface conditions, as well as human activities (Fig. 1c; 
Wang et al., 2021). Generally, the concentrations of hydrochemical components (TDS, CI, S04”, 
Na’, K’, Ca”, and Mg”) in the shallow groundwater within the oasis showed horizontal zoning 
and increased along the regional flow direction from northwest to southeast (Fig. 6). The TDS of 
shallow groundwater samples ranged from 1021.9 to 24,134.6 mg/L. The high value of dissolved 
ion concentrations for shallow groundwater was observed in several sporadic areas in the 
southeastern and central parts, whereas the low value was distributed in the northwestern part. 
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Overall, the spatial distribution pattern of dissolved ion concentrations in shallow groundwater 
was consistent across the oasis but lower in the northwest than in the southeast (Fig. 6), which 
could be explained by the groundwater depth, dissolving mineral, agricultural irrigation, and 
drainage in this region (Pant et al., 2018). 
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Fig. 6 Spatial distributions of shallow groundwater hydrochemistry in the Wei-Ku Oasis. (a), TDS; (b), CI; (c), 
SO,” (d), Na"+K*; (e), Ca”’; (P), Mg”. 
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Fig. 7 Spatial distributions of middle groundwater hydrochemistry in the Wei-Ku Oasis. (a), TDS; (b), CI; (c), 
S04”; (d), Na*™+K"; (e), Ca”; (f), Mg”. 
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Furthermore, as shown in Figure 7, the spatial variation trend of middle groundwater 
hydrochemistry was not obvious in the oasis, and the spatial distribution pattern differed between 
various ions, probably due to the spatial heterogeneity of lithology and groundwater exploitation 
(Fig. 1c) (Wang et al., 2013). The TDS value of middle groundwater ranged from 561.0 to 2934.8 
mg/L. The spatial distribution patterns of TDS, CI, and SO,” of middle groundwater were 
consistent and were much higher in the central part than in the surrounding areas (Fig. 7). 
However, the cation concentrations (Na’, K*, Mg’, and Ca’) of middle groundwater did not 
exhibit an obvious spatial variation trend. In addition, for middle groundwater within the oasis, 
the recording sites with the highest and/or lowest value varied significantly among various ions. 
Among the major cations, the Ca** concentration of middle groundwater was higher in the 
western part of the oasis than in the eastern part, with the maximum value observed near the 
Xinhe County. Compared to other regions, the Mg” concentration was much higher in the central 
and western parts while the Na’ and K“ concentrations in the central and southwestern parts were 
found relatively high. 


3.3 Seasonal variations of hydrochemistry 


Figure 8 exhibited the seasonal variations of the hydrochemical components (TDS, CT, SO,’ , 
Na’, K“, Ca”*, and Mg”) for surface water and groundwater in our study region. The dissolved 
ion concentrations of surface water and groundwater fluctuated during the sampling period, 
showing significant seasonal variations (Fig. 8). Usually, except for K”, the seasonal variation 
patterns of major dissolved ion concentrations (Ca”, Mg”, Na’, Cl, and SO,” ) for river water 
samples in the upstream were consistent. The lowest values of Cl, S04”, Nat, Mg”, and Ca” 
concentrations of river water were observed in summer (26.7, 75.5, 23.9, 15.7, and 72.2 mg/L, 
respectively), whereas the highest values were detected in spring (67.4, 162.7, 42.2, 25.0, and 
87.0 mg/L, respectively). Nonetheless, the highest and lowest values of K* concentration were 
observed in autumn and spring, respectively (Fig. 8e). In contrast, the seasonal differences of CI , 
SO,’ , Na’, and Mg” concentrations of river water were significant, while the values of K* and 
Ca™ concentrations varied slightly among seasons. 

As shown in Figure 8, the ion concentrations of middle groundwater were much lower than 
those of shallow groundwater across the Wei-Ku Oasis in each season. Hydrochemical 
components of shallow and middle groundwater samples varied significantly in different seasons, 
revealing the seasonal variations of hydrochemical components, which may be related to the 
water sources and water-salt interaction between different aquifers (Zeng et al., 2020). Except for 
Na’ and Mg”, the seasonal variation patterns of major dissolved ions (Cl , S04”, K*, and Ca”) 
of middle groundwater in the oasis were consistent. However, the seasonal variations of different 
ion concentrations of shallow groundwater differed (Fig. 8). For the middle groundwater, the 
values of Cl and Ca?” concentrations were slightly higher in spring (CI=280.9 mg/L and 
Ca**=160.4 mg/L) than in summer (Cl =258.7 mg/L and Ca**=140.4 mg/L), whereas the Na’ and 
Mg” concentrations were lower in spring (Na™=175.6 mg/L and Mg”’=80.7 mg/L) than in 
summer (Na’=232.6 mg/L and Mg?*=88.2 mg/L). For shallow groundwater, the highest and 
lowest values of Cl, SO,”, and Na’ concentrations were observed in spring and summer, 
respectively. Meanwhile, the K* and Ca” concentrations were highest in autumn (Ca”=5 14.1 mg/L 
and K*=28.9 mg/L) and lowest in summer (Ca”"=397.5 mg/L and K*=18.1 mg/L), and the Mg” 
concentration was highest in spring (358.3 mg/L) but lowest in autumn (309.0 mg/L) (Fig. 8). 


4 Discussion 


4.1 Spatiotemporal variation characteristics of hydrochemistry 


The mean ion concentrations (Ca”, Mg”, Cl, and SO,” ) of river water in our study area were 
much higher than the global average (Meybeck, 2003), but lower than those in the Tarim River 
(Table 2), probably reflecting the lower precipitation and stronger evaporation in the Tarim River 
Basin (Pant et al., 2018). Similarly, the Ca”, Cl, and SO, concentrations of river water in the 
Wei-Ku Oasis were higher than those in the Shule River, Heihe River, and Shiyang River, while 
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Fig. 8 Seasonal variations of surface water and groundwater hydrochemistry in spring, summer, and autumn 
across the study area. (a), TDS; (b), CI; (c), S04”; (d), Na’; (e), K*; (f), Ca™*; (g), Mg’. IQR, interquartile range. 


Mg” concentration of river water in the Wei-Ku Oasis was higher than that in the Shule River but 
lower than those in the Heihe River and Shiyang River (Table 2). Such findings may be attributed 
to the different water sources and evaporation (Wang et al., 2013; Liu et al., 2018). The 
concentrations of dissolved ions of groundwater samples displayed wide ranges, indicating that 
groundwater had undergone different hydrogeochemical processes within the aquifer system (Liu 
et al., 2018). Another reason may be the water exchange between confined water and higher 
salinity phreatic water in the oasis (Zeng et al., 2020). In recent decades, farmland expansion and 
population increase have resulted in a sharp increase in water consumption. Due to the limited 
surface water resources, water consumption for agricultural irrigation was increased and 
groundwater level was decreased, which in turn caused phreatic water leakage downward into the 
confined aquifer in response (Wang et al., 2021). Moreover, the ion concentrations of water 
samples significantly varied among different water components and were much higher in shallow 
groundwater than in river water and deep groundwater (Table 1). The higher concentrations of 
dissolved ions in shallow groundwater samples may be attributed to the fact that high temperature 
and surface water infiltration could promote the erosion, dissolution, and leaching processes of 
soluble minerals within the soil profile, causing large amounts of soluble salts to migrate to the 
shallow aquifers (Tweed et al., 2011; Ma et al., 2013; Huang et al., 2017). On the other hand, a 
large amount of liquid and solid wastes from nearby cropland, settlements, and factories were 
directly or indirectly discharged into the aquifers, which critically affected the hydrochemical 
components of groundwater in the oasis (Sharma et al., 2012; Pant et al., 2018). 
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Table 2 Hydrochemical compositions of river water for the Weigan-Kuqa River Basin and their comparisons 
with other rivers in arid areas of China and the global mean 


River cs Me pis ae ae SOE Reference 
mg/L 
Weigan-Kuqa River 79.2 20.7 37.4 145.9 46.2 122.9 This study 
Shule River (China) 28.8 19.3 13:9 139.5 11.1 34.8 Zhou (2015) 
Heihe River (China) 57.6 21.5 18.1 204.4 10.0 86.1 Nie et al. (2005) 
Shiyang River (China) 65.7 29.3 4.9 64.9 1.8 30.9 Gao et al. (2006) 
Tarim River (China) 56.7 145.3 302.4 - 769.3 295.1 Wang et al. (2013) 
Global mean 15.0 4.1 8.6 58.4 7.8 11.2 Meybeck (2003) 


Note: -, no data. 


Spatial heterogeneity of groundwater hydrochemical components was evident in the oasis 
(Figs. 6 and 7), which probably overlapped with soluble minerals, groundwater depth, and 
cropland (Figs. 1 and 9) (Pant et al., 2018). Most major ion concentrations of shallow 
groundwater within the Wei-Ku Oasis followed the same increasing trend along the flow direction 
from northwest to southeast (Fig. 6). The higher concentrations of shallow groundwater 
hydrochemistry in the southeast of the Wei-Ku Oasis were found mainly in the desert region with 
relatively low groundwater depth (Fig. 9); it is the area most affected by strong evaporation 
(Huang et al., 2017). In contrast, lower concentrations of shallow groundwater hydrochemistry in 
the northwestern oasis were considered to be the outcomes of irrigation and drainage processes as 
well as groundwater extraction (Thomas et al., 2015), due to stronger dilution by agricultural 
irrigation water and more salt discharge by agricultural drainage canals (Wang et al., 2018; Wang 
et al., 2021). Moreover, a higher TDS value of middle groundwater was found in the central part 
of the Wei-Ku Oasis (Fig. 7), which may be related to land use types and groundwater 
overexploitation. Cultivated area expansion has caused groundwater overexploitation in the arid 
oases (Huang et al., 2017; Wang et al., 2021). Groundwater overexploitation could alter the 
leaking direction between aquifers and redistribute soluble salt within the aquifers (Zeng et al., 
2020; Wang et al., 2021). 
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Fig. 9 Spatial distributions of land use and land cover (LULC; a) and groundwater level depth (b) in the Wei-Ku 
Oasis in 2018. The data of groundwater level depth were from Wang et al. (2021). 


Concentrations of most major ions in groundwater varied significantly among each season in 
the arid oasis (Fig. 8); they were closely related to air temperature, runoff, groundwater depth, 
and agricultural activities (Wang et al., 2013), suggesting that the seasonal variation 
characteristics of groundwater salt concentration can reflect the temporal variations of recharge 
sources and methods, as well as hydrological and hydrogeochemical processes (Beal et al., 2019; 
Wang et al., 2021). The hydrochemical values of river water in the upstream were lower in 
summer and higher in spring, which contradicts the seasonal variations of rainfall and runoff. 
Furthermore, the heavy runoff (Fig. 2) generated from the coupled effects of intense glacier/snow 
melting and high rainfall during the wet season was responsible for the mild low concentrations 
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of hydrochemical components in surface water (stronger dilution effect) (Zhou, 2015; Liu et al., 
2018). Additionally, the TDS value of shallow groundwater was higher in spring than in summer 
and autumn, which may be ascribed to the infiltration from agricultural irrigation water, following 
the presence of chemical fertilizers and pesticides, the occurrence of dissolving mineral and soil 
salinity, and the dumping of post-harvest residues (flood irrigation in winter and spring) (Sharma 
et al., 2012; Wang et al., 2021). Furthermore, a higher TDS value of middle groundwater was 
observed in summer (Fig. 8), which was closely related to the cropland area and groundwater 
level depth in summer (Fig. 9), suggesting the accumulated influences of human activities 
(agricultural irrigation and groundwater extraction) (Pant et al., 2018; Castellano et al., 2019). 


4.2 Dominant factors of groundwater hydrochemistry evolution 


The diagram of the weight ratios between TDS and Na‘/(Na‘+Ca’’) or Cl /(CI+HCO; ) exhibited 
the major mechanisms controlling the hydrochemistry of nature water bodies (Gibbs, 1970), 
including rock dominance, precipitation dominance, and evaporation-crystallization process 
(Wang et al., 2013) (Fig. 10). As shown in Figure 10, the weight ratio Na’ /(Na‘+Ca’’) of river 
water in the upstream of the Wei-Ku Oasis was lower than 0.5, and the salinity was low, 
indicating that the hydrochemical compositions of river water is primarily controlled by rock 
dominance. Also, referring to natural mechanisms controlling the groundwater hydrochemistry, 
there was an obvious difference between various aquifers in the Wei-Ku Oasis. It was noted that 
the hydrochemical compositions of phreatic water were mainly controlled by rock dominance and 
evaporation-crystallization process, likewise to the deep groundwater hydrochemistry that was 
primarily influenced by rock dominance mechanism. Still, the weight ratio Na‘/(Na‘+Ca’’) of 
some middle groundwater samples was higher than 0.6 (Fig. 10), suggesting that additional 
mechanisms contribute to the hydrochemical compositions of middle groundwater in the arid 
oasis (e.g., human activities) (Zeng et al., 2020). The major cations and TDS of deep groundwater 
in the oasis were similar to those of confined groundwater in the upstream of the oasis, which 
again reflected the deep groundwater hydrochemistry to be majorly controlled by rock dominance 
mechanism (such as interactions between water and rock) (Ma et al., 2013; Pant et al., 2018). 
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Fig. 10 Plots of TDS vs. Na‘/(Na‘+Ca’") in surface water (a, b) and TDS vs. CIl’/(CI+HCO; ) in groundwater 
(c, d) in the Wei-Ku Oasis. The basis for this plot was from Gibbs (1970). 


Furthermore, dissolved species and their relations could reveal the solute origin and 
hydrogeochemical processes (Fisher and Mullican, 1997; Wang et al., 2013). As presented in 
Figure 4, river water samples were mainly located near the Ca end-member and the HCO3+CO3 
end-member, informing that the major ionic compositions were mainly derived from carbonate 
dissolution (Xiao et al., 2015). In general, the hydrochemical compositions of deep groundwater 


JOURNAL OF ARID LAND 


samples were mainly derived from carbonate dissolution and evaporite, meanwhile the 
hydrochemical compositions of middle groundwater samples came from the carbonate 
dissolution, and those of shallow groundwater were from the evaporite dissolution in the Wei-Ku 
Oasis (Pant et al., 2018; Zeng et al., 2020). Additionally, several geochemical processes such as 
evaporation-condensation, dissolution and precipitation, cation exchange/adsorption, and mixing 
would be responsible for different groundwater types in the Wei-Ku Oasis (Wang et al., 2013; Li 
et al., 2019). Phreatic groundwater salinization in the arid oasis may be resulted from the 
increasing ionic concentrations due to the evaporation-condensation of recharge water (infiltration 
from riverbank, canal, reservoir, and irrigation), as well as the interactions between rock and 
water. This eventually leads to the mineral deposition, forming saturated zones (Wang et al., 
2021). Furthermore, the Na‘/Ca”* ratio of groundwater samples in the oasis increased along the 
flow path, implying the Ca” loss due to carbonate precipitation and the additional Na‘ increase 
through cation exchange (Tweed et al., 2011; Liu et al., 2018). As well, the mixing process 
between groundwater and environmental flows within the aquifer system could also play an 
important role in groundwater hydrochemistry, especially for shallow and middle groundwater, 
including the mixing between groundwater and weakly mineralized surface water, as well as 
between confined water and high mineralized phreatic water in some areas (Su et al., 2009; Zeng 
et al., 2020). 

Human activities such as land use changes, groundwater extraction, agricultural irrigation, and 
drainage have significantly affected the hydrochemistry of phreatic water and confined water in 
the Wei-Ku Oasis in recent decades (Xiao et al., 2015; Li et al., 2019). In the Tarim River Basin, 
evapotranspiration results in the accumulation of evaporites in the upper layer of the soil profile, 
and the inappropriate irrigation regime can exacerbate soil salinization, thus leading to the 
widespread soil salinization in our study region (Liu et al., 2018; Chen et al., 2019). In parallel, 
Wang et al. (2021) reported that soil salinization in the irrigated district was notable due to intense 
evaporation, long-term irrigation, and relatively shallow groundwater in the Wei-Ku Oasis; they 
also observed that the soil salt content in the 0-100 cm soil depth was 3.98 g/kg in 2017. Our 
findings revealed that the cropland area in the oasis is depended on river water and groundwater 
extraction for irrigation purposes (Table 3) and has expanded over decades (Fig. 9). Consequently, 
these activities resulted in the irrigation return flow as a dominant recharge source of shallow 
groundwater. 


Table 3 Information related to land cover types, human activities, and groundwater level in the Wei-Ku Oasis 
from 2000 to 2015 


Parameter 2000 2010 2015 
Cropland area (km?) 3969.3 4302.8 4389.5 
Bare land area (km’) 1432.7 1077.8 965.3 
Water-saving irrigation area (km’) 0.0 460.5 1103.0 
Irrigation water amount (x10* m°) 22.3 25.3 23.1 
Salt recharged by inflow water (*10* t) 89.0 101.3 92.6 
Agricultural water discharged by drainage canals (x10° m°) 3.1 23 1.9 
Agricultural salt discharged by drainage canals (x 10* t) 170.6 150.9 94.5 
Groundwater exploitation (x10° m°) 0.5 2.8 2.9 
Depth to groundwater level (m) 2.7 4.0 5.7 


Note: The data of human activities and groundwater level were from Wang et al. (2021). 


The vertical infiltration from irrigation water could promote the dissolution and leaching of soil 
salt and soluble mineral by agricultural activities, thus transferring the soluble salt into shallow 
groundwater and increasing ion concentrations (Han et al., 2011). Previous studies observed that 
groundwater salinity had a significant linear correlation with soil salinity in agricultural and 
saline-alkaline land in arid inland oasis areas (Wang et al., 2021), indicating that irrigated 
cultivation in arid inland oasis areas could easily affect groundwater salinity owing to intensive 
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irrigation and inefficient drainage. Likewise, the improved agricultural drainage channels could 
also significantly affect groundwater hydrochemistry, and shallow groundwater with high salinity 
near cultivated land will be discharged through the drainage canals (Table 3), thus decreasing the 
groundwater level and salinity (Wang et al., 2021). A large amount of salt was discharged by 
drainage canals in the irrigation district of the Wei-Ku Oasis, with a value of 94.5x10* t in 2015 
(Table 3). Furthermore, due to river water deficiency, groundwater (especially confined water) 
was overexploited for agricultural irrigation during the growing season over decades, causing a 
reduction in groundwater level (Table 3). Therefore, groundwater exploitation intensifies the 
process of groundwater salinization in this oasis to some extent (Wang et al., 2013), leading to the 
downward leakage of phreatic groundwater into confined water with the decline of groundwater 
level. The process also promotes the mixing and interaction of fresh water and saline water 
among various aquifers (Zeng et al., 2020; Wang et al., 2021). Relevant researches have shown 
that shallow groundwater depth and intense evaporation are the primary factors influencing the 
spatial distribution of phreatic groundwater hydrochemistry (Wang et al., 2021). 


4.3 Implications for groundwater salt migration under the influence of human activities 


As mentioned above, spatiotemporal distributions of groundwater hydrochemistry within the 
aquifer system in arid inland oasis areas would provide insight into the source of groundwater salt 
and its migration mechanism under the influence of human activities (Wang et al., 2021). 
Groundwater salt migration is influenced by many factors in the arid dried-up river oasis, such as 
climate conditions (precipitation and temperature), runoff, hydrogeological conditions, land cover 
types, groundwater extraction, agricultural irrigation, and drainage (Pant et al., 2018; Porhemmat 
et al., 2018; Castellano et al., 2019). In the upstream of the Weigan-Kuga River Basin, salt ions in 
river water migrate to the downstream with the migration primarily controlled by rock 
dominance. In the oasis, river water flows into the field for agricultural irrigation through canals 
and experiences strong evapotranspiration before infiltration, but salt ions in river water can also 
transfer downward into the soil layer (Chen et al., 2019). Within the vadose zone, soil moisture 
supplied by the infiltration of rain, riverbank, canal, reservoir, and irrigation could migrate 
downward to recharge phreatic groundwater by gravity gradient, while soluble ions in soil 
moisture also migrate downward and accumulate in the phreatic aquifer due to evaporation, 
surface water dilution, soil salt dissolution, leaching, and cation exchange (Xiao et al., 2015). Soil 
moisture and salt patterns in the irrigated district were reshaped by long-term irrigated agriculture 
(Yin et al., 2021). Within the aquifer, phreatic groundwater supplied by surface water infiltration 
and lateral groundwater flow could leak downward to confined aquifer under the driving of water 
potential gradient due to groundwater overexploitation. Meanwhile, high concentrations of 
dissolved ions in phreatic groundwater also transfer downward into the confined aquifer, leading 
to the deterioration of confined water quality due to the interaction of saline water and fresh water 
that are mainly controlled by the evaporation-crystallization process, rock dominance (strong 
water-rock interaction), and mixing mechanism. So, large-scale land reclamation and groundwater 
extraction have dramatically altered groundwater water-salt migration processes in this region. 


5 Conclusions 


Spatiotemporal variations and evolution of groundwater hydrochemistry were investigated using 
data of 199 water samples collected in the Wei-Ku Oasis of Tarim Basin. The results showed that 
the mean ion concentrations of river water in the Wei-Ku Oasis were much larger than the global 
average but lower than the Tarim River. Groundwater hydrochemistry showed spatiotemporal 
dynamics and seasonal differentiation, and its spatial distribution was complex. Vertically, 
groundwater ion concentrations decreased with depth along with the aquifers. Generally, the TDS 
and 8'*O values of river water in the upstream increased from west to east, while shallow 
groundwater ion concentrations in the oasis increased from northwest to southeast. Higher values 
of TDS were found in spring for shallow groundwater while in summer for middle groundwater. 
Pronounced spatiotemporal heterogeneity demonstrated the impacts of geogenic, climatic, and 
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anthropogenic conditions. Agricultural irrigation and drainage, land cover change, and 
groundwater extraction reshaped the spatiotemporal patterns of groundwater hydrochemistry. 
Lately, groundwater overexploitation altered the leaking direction between the phreatic and 
confined aquifers, causing the salt-fresh water interaction and groundwater environment 
deterioration. 
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